Magnetic anisotropy and reversal in epitaxial Fe/MgO(001) films revisited 
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We investigate the magnetization reversal in Fe/MgO(001) films with fourfold in-plane magnetic 
anisotropy and an additional uniaxial anisotropy whose orientation and strength are tuned using 
different growth geometries and post growth treatments. The previously adopted mechanism of 180° 
domain wall nucleation clearly fails to explain the observed 180° magnetization reversal. A new 
reversal mechanism with two successive domain wall nucleations consistently predicts the switching 
fields for all field orientations. Our results are relevant for a correct interpretation of magnetization 
reversal in many other epitaxial metallic and semiconducting thin films. 

PACS numbers: 75.60.Jk, 75.30.Gw, 75.70.Ak 



Magnetic anisotropy is one of the most important 
properties of metallic and semiconducting thin-film mag- 
nets [H, 0] . In magnetic films of cubic systems an in- 
plane fourfold magnetic anisotropy is expected, but of- 
ten an additional uniaxial magnetic anisotropy (UMA) 
is observed to be superimposed on top of the fourfold 
anisotropy 0, 0| • The extra UMA has been attributed 
to different origins, including a self-shadowing effect oc- 
curring during oblique deposition [5j, [y, [7( , the bonding 



mm 



and the Neel surface 
ll| . Moreover, ion sput- 



between film and substrate 
effect on a stepped substrate 
tcring has been demonstrated as a reliable tool to control 
the orientation and streng th of UMA [H [H. 

When applying thin films for magnetic data storage 
and spintronic devices, the magnetization reversal mech- 
anisms and their dependence on the anisotropy symme- 
try need to be known and controlled. The magnetiza- 
tion reversal process for combined cubic and uniaxial 
anisotropies is sensitive to the specific anisotropy geom- 
etry and strength [hi . [TBI . Depending on the field ori- 
entation, hysteresis curves with one and two steps are 
observed in various films, and explained in terms of nu- 
cleation and propagation of 90° and 180° domain walls 
(DWs) [H, E3]. A model based on minimizing the mag- 
netic energies has been introduced with DW nucleation 
energies ego° for 90° DWs and eiso° for 180° DWs, re- 
spectively, in order to account for the observed switching 
fields [l7| . A special magnetic switching process involv- 
ing three steps can be observed when the additional UMA 
along the cubic easy axis exceeds ego<> d8, 19]. Until now, 
such a switching has been assumed to be mediated by two 
90° DW nucleations at the first and the third step and 
one 180° DW nucleation occurring in between. 

Previous work on Fe/MgO(001) films grown at normal 
incidence revealed a weak UMA along Fe(010) [joj . Re- 
cently, we successfully relied on ion sputtering to manip- 
ulate the strength of the in-plane UMA along Fe(110) 
in Fe/MgO(001) [HI]. Pai ' k et al. found that a pro- 
nounced UMA can be induced in Fe/MgO(001) by rely- 



ing on oblique-incidence molecular beam epitaxy (MBE) 
growth [221] . Up to now, all of the measured hysteresis 
loops in Fe/MgO(001) only revealed one or two steps. 

Here, we report on a detailed study of magnetization 
reversal in Fe/MgO(001) films, where strength and ori- 
entation of UMA are tuned either by ion sputtering or by 
oblique- incidence MBE growth. A novel mechanism is in- 
troduced with two successive DW nucleations to explain 
the 180° magnetic switching that occurs for one-step and 
three-step loops. Our model consistently explains the ex- 
perimental results for films with different UMA, revealing 
the universal nature of the magnetization reversal. 

In general, the in-plane UMA, which is superimposed 
on the cubic anisotropy K\ of Fe, can be separated 
into two components: K u \ along [010] and K U 2 along 
[110] [23]. If K u i <C Ki and K u2 < K 1} the com- 
ponent K U 2 rotates the position of the overall easy 
axes backwards with respect to the uniaxial hard axis 
over an angle 5 that is approximately given by S = 
■i sin^ 1 (K u2 / Ki) [HI' as illustrated in Fig. 1(a). In order 
to obtain a UMA with different orientation and strength, 
three Fe films were grown by MBE on MgO(OOl) sub- 
strates using different growth conditions and post growth 
treatments. For sample A, the incident Fe beam was at 
an angle of 49° with respect to the surface normal and 
with azimuthal angle along Fe[010]. During deposition 
of sample B at normal incidence, the substrate was ro- 
tated around the surface normal. The nominal thickness 
of samples A and B was 15 nm, as monitored by a cal- 
ibrated quartz crystal oscillator. The growth geometry 
for sample C was the same as for sample B, but the nom- 
inal thickness was 100 nm. Subsequently, sample C was 
sputtered with 2 keV Ar + ions at an incidence angle of 
60° with respect to the surface normal and with azimuth 

fixed in between Fe[100] and Fe[l 10] . After sputtering for 
250 minutes, the film thickness was reduced to 15 nm, as 
verified by ex situ x-ray reflectometry. Before removing 
the samples from the vacuum chamber, they were capped 
with a 2 nm thick protective Au layer. The magnetic 
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properties were measured by the longitudinal and trans- 
verse magneto-optical Kerr effect (MOKE) for different 
field orientation <f> as denned in Fig. 1(a). 

For sample A, a considerable UMA along [010] is intro- 
duced by the oblique deposition. Three-step loops as well 
as one-step and two-step loops are observed at different cj>, 
as illustrated in Figs. 1(b) to 1(d). The switching events, 
which occur for increasing field and 0° < <j> < 90°, are 

[100]->[100] for the one-step loops, [010]->[100]->[010] 

for the two-step loops, and [010]->[100]-v[100]->[010] for 
the three-step loops [19j . The magnetization switches by 
180° for the one-step loops and for the second step of 
the three-step loops, and by 90° for the other steps. The 
corresponding spin orientations are marked by the arrows 
that are enclosed in a square in Figs. 1(b) to 1(d). 

Up to now 90° as well as 180° DW nucleations have 
been invoked to interpret the 90° and 180° magnetic 
transitions, respectively [It], El- The coercivity re- 
lated to the DW nucleation energy can be derived from 
the energy gain between the local minima at the ini- 
tial and final easy axes involved in the transition 17 1. 
The theoretical switching fields are obtained as H c \ = 
(cqqo — K u x)/[M(smcf) — cos0)] for the magnetic switch- 
ing process [010]-»[100], H c2 = (e 90 ° + K ul )/[M(sm.<f> + 

cos0)] for [100M010], H c3 = (e 90 o + K ul )/[M(smcf> - 
cos</>)] for [100M010], H ci = (e 90 o - K ul )/[M(cos - 
sin</>)] for [010]->[100], and H c = ei 8O o/[|2M(cos0)|] for 

[100]— >[100], where M is the magnetization. 

For sample A, the <fi dependence of the measured 
switching fields can be nicely fitted using the theoreti- 
cal switching fields (see Fig. 2(a)), provided we assume 
the switching fields correspond to H c i, H c2 and H c3 (see 
Figs. 1(c) and 1(d)). The fitting results in the parameters 
K ul /M = 2.70 ± 0.02 mT and e 90 o/M = 0.61 ± 0.02 mT, 
where K u \ > e 9 o° is the necessary condition for the occur- 
rence of three-step loops 



17. 1 



wo 



[19J. Following Refs. 
try to describe the experimental switching field iJ c igo°, 
which corresponds to a 180° magnetic transition, in terms 
of 180° DW nucleation. Using the corresponding theo- 
retical switching field H c , we obtain the fit that is shown 
in the two insets to Fig. 2(a). H c reaches a minimum 
at <f> = 0° an d a maximum at (j> = 90°. The experi- 
mental switching field -ff c i8o° reveals, however, a peak at 
both <j> = 0° and 90°. Moreover, around <j> = 90° H c has 
a slope that disagrees with experiment. 180° DW nucle- 
ation clearly fails to describe the 180° magnetic transition 
in sample A. Surprisingly, the theoretical expression for 
H c2 , which corresponds to a 90° DW nucleation, allows 
us to fit the -ffci80° data (see below for more details). 

The universal character of the absence of 180° DW 
nucleation for Fe/MgO(001) is confirmed by the <f> de- 
pendence of the switching fields for samples B and C in 
Figs. 2(b) and 2(c). Only one-step and two-step loops are 



switching route for the two-step loops appearing in sam- 
ples B and C (180° -<J -> 90°+<5 -> -6 for -6 < (/>< 45°) 
is different from the path for sample A. When 6 = 0, i.e., 
K u2 = 0, the experimental switching fields for this type 
of two-step loop should correspond to the above derived 
expressions for H c2 and H C 4. In case 6^0, i.e., K u2 ^ 0, 
the theoretical expressions need to be extended: 



H c2 = 



e 9 o°-25 + if «i (cos 2 S - sin 2 6) 
M[cos(0 + 6) + sm(0 - 6)} 

e 9 o°+2<5 - K ul (cos 2 6 - sin 2 6) 
M[cos(0 + 6)- sm((j) - 6)] 



observed. As discussed in our recent publication [21(, the 



where e 9 o°-2i5 and e 9 o°+2i5 are the corresponding DW nu- 
cleation energies. 

For sample B the two switching fields (H c2 , H c ^) 
have a dependence on <p that is symmetric about (100). 
Moreover, the angular dependence of H c2 reveals a 
clear, abrupt step when crossing (110), as illustrated in 
Fig. 2(b). We conclude that sample B has a small in- 
plane UMA along [010] [17| . The fitting parameters are 
K ul /M = 0.19 ± 0.01 mT and e 90 =/Af = 0.36 ± 0.01 mT. 
Because K u \ < e 9 o° , three-step loops cannot be observed. 

In sample C, UMA with components along both [010] 
and [110] is introduced by the Ar + ion sputtering. The 
overall easy axes are observed to deviate from (100) by 
an angle 8 = 3°, i.e., K u2 jK\ « 0.1. From the results in 
Fig. 2(c) we find that the UMA component along [010] 
and the DW nucleation energies are K u \jM = 1.69 ± 
0.02mT, e 90 o_ 2S /M = 1.83 ± 0.02 mT, and e 90 ° +2S /M = 
2.29 ± 0.02 mT, respectively. Because K u \ is comparable 
to e 90 °-2,5, the path 270° + S -> -5 -> 90° + S is ener- 
getically more favorable when compared to the counter- 
clockwise path via 180° — S for the whole range of angles 
45° < 4> < 135°. Consequently, both H c2 and H C 4 change 
monotonously within this range. 

We again try to describe -ff c i80° for the one-step loops 
of samples B and C in terms of 180° DW nucleation, as il- 
lustrated in the insets to Figs. 2(b) and 2(c), respectively. 
The theoretical curves clearly disagree with experiment. 
From the fitting for our samples with different anisotropy 
geometry and strength we conclude that 90° ± 25 mag- 
netic transitions in Fe/MgO(001) films are mediated by 
90° ± 26 DW nucleation, but 180° magnetization reori- 
entations are not mediated by 180° DW nucleation. 

Which mechanism dominates the 180° magnetic rever- 
sal? Based on the obtained values for K ul and e 9 o°, we 
plot the energy difference between the relevant easy axes 
as a function of the applied field in Fig. 3(a) for sample 
B at 4> = 10° and in Fig. 3(b) for sample A at (j> = 65°, 
respectively. In the previously adopted model, the en- 
ergy difference between 180° and 0° is treated in terms 
of one single barrier, and 180° DW nucleation occurs 
at H c when A_Eigo°^o° = e i80°, where eiso° is assumed 
to correspond to 2e 9 o° [19]- According to our analysis 
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the switching between [100] and [100] is governed by two 

separate energy barriers between [100] and [010] and be- 
tween [010] and [100], respectively. The switching then 
corresponds to two 90° DW nucleation processes. The 

energy barrier for the transition from [100] and [010] 
becomes Ai?iso°^9o° = e 90° a t -ff C 2- However, since 
Ai59o°_»o o already exceeds ego° at H C 2, the domains along 
[010] are unstable and cannot grow. Therefore, a second 
nucleation of domains along the final [100] remanent di- 
rection occurs at H c2 , and the two successive 90° DW nu- 
cleations appear as one single step in the MOKE loops. 
In case 5^0, this process consists of a 90° — 28 DW 
nucleation and a subsequent 90° + 25 DW nucleation, or 
vice versa. Based on our new model the experimental 
switching fields -ff c iso° for all three samples are fitted by 
the expressions for H C 2 in Figs. 2(a) to 2(c) and the in- 
sets. This way, all switching fields can be nicely fitted by 
consistently using the same £9o°±2<5 and K u \ values for 
the complete range of angles. 

In case of two successive DW nucleations, H c i is not 
an experimental observable switching field. H C 4 only 
indicates AEqqo + s^_s = 6900+25. We have plotted in 
Figs. 2(a) and 2(b) the "virtual" H C 4 values for samples 
A and B. When 0° < (j> < 45°, the two successive DW 
nucleations appear as one-step loops for H C 2 > H C 4 (see 
Fig. 3(a)). For H C 2 < H C 4, the magnetization loops re- 
veal a two-step behavior with two separate 90° DW nucle- 
ations at H C 2 and H C 4, respectively. When 45° < <j) < 90° 

and K u i > £90°, the magnetization switches from [010] 

to [100] at H cl , where A_E 2 7o°^i80° = e 90°- Because 
AEg a^ a decreases with increasing applied field and be- 
comes ego° at H C 4 (see Fig. 3(b)), two successive DW 
nucleations appear for H c2 < H C 4- Upon further increas- 
ing the field, A_Ego°_»o° becomes negative, and finally 
reaches —£90° at H c ^, where the magnetization switches 
backwards from [100] to [010]. As a result, the magne- 
tization loops contain three steps. For H C 2 > H C 4, the 
domains aligned along 90° are energetically stable when 
the applied field exceeds H c2 , resulting in two-step loops. 

By comparing the expressions for H C 2 and H C 4 the field 
orientation for the occurrence of one-step or three-step 
loops can be obtained. For the simple case 5 = 0, the 
condition tan</> < if M i/ego° needs to be satisfied, where 
< 4> < 45° for the one-step loops and 45° < cp < 90° 
for the three-step loops, respectively. Our model predicts 
that the ranges of angles for which a one-step loop should 
be observed are -45° < < 45°, -28° < < 28°, and 
—43° < <\> < 36° for samples A, B, and C, respectively. 
The critical angles separating the occurrence of two-step 
and three-step loops are 4* = 90° ± 13° for sample A. Our 
model calculations nicely agree with experiment. 

Introducing two successive DW nucleations allows us 
to consistently interpret the 180° magnetization reorien- 
tation. Since the DWs induced during the first nucleation 



are energetically unstable, it should be very hard to ob- 
serve these intermediate domains. In real films, however, 
DW nucleation and propagation are often perturbed by 
defects and roughness [24(, implying magnetic switching 
is not as sharp as predicted. As illustrated in Fig. 4(a), 
we observed the metastable intermediate states in the 
loops of sample A close to the critical angles separating 
the occurrence of two-step and three-step loops. The blue 
curve (increasing field) reveals an overshoot for magnetic 

switching from [100] to [100], indicating that the Fe spins 
align for a short time along [010] before jumping to [100]. 
The red curve (decreasing field) reveals a similar feature. 
We note that the second intermediate state in the blue 
curve in Fig. 4(a) is not collinear with the first intermedi- 
ate state in the red curve and vice versa, which is different 
from the loop shown in Fig. 1(d). The non-collinearity 

implies that not all spins switch from [100] to [010] and 
then to [100], but some of the spins remain aligned along 
[010]. This points to coexistence of magnetic switching 
processes with both two-step and three-step loops. The 
overshoot can still be observed a few degrees away from 
the critical angles. In sample B we observe a mixture of 
one-step and two-step loops, as illustrated in Fig. 4(b). 
The red curve of the transverse MOKE loop reveals two 
separate 90° DW nucleations, while the blue curve cor- 
responds according to our model to two successive and 
indistinguishable DW nucleations. Experimentally, the 
transition between two reversal mechanisms does not oc- 
cur at one critical angle as predicted by theory, but ex- 
tends over a small finite range of angles. Time resolved 
MOKE may be able to detect the ultrafast magnetization 
dynamics and the intermediate domain formation in the 
process of two successive DW nucleations. 
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FIGURE CAPTIONS 

FIG. 1. (Color online) (a) Definition of the angles that 
are used to describe a film with in-plane cubic anisotropy 



and UMA. Typical longitudinal MOKE loops for sample 
A with (b) one step at (f> = 8°, (c) two steps at 4> = 88° 
and (d) three steps at (f) = 68°. The blue (red) curves 
are for applied fields varying from negative (positive)to 
positive (negative) saturation. The arrows enclosed by a 
square represent the orientation of the Fe spins. 

FIG. 2. (Color online) The experimental switching 
fields (symbols) as a function of the field orientation (j>, 
and the corresponding theoretical curves for H c \ (ma- 
genta), H c2 (red), H c3 (purple), H C 4 (blue), the "virtual" 
H C 4 (dashed gray), and H c (green) for (a) sample A, (b) 
sample B, and (c) sample C. 

FIG. 3. (Color online) Energy differences AEi SO °^o° 
(red), A_E 2 70°— iso° (cyan), A£;i 80 »^9 o (green), and 
A-Ego°_^o° (blue) as a function of the applied field for 
(a) sample B at <j> = 10° and (b) sample A at cf> = 65°. 

FIG. 4. (Color online) (a) The longitudinal magnetiza- 
tion loop of sample A at <j> = 76° and (b) the transverse 
magnetization loop of sample B at <f) — 29° reveal over- 
shoots near the critical angles for the field orientation. 
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